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Abstract. The luminous infrared-loud quasar IRAS 13349+2438 was observed with the XMM-Newton Observatory 
as part of the Performance Verification program. The spectrum obtained by the Reflection Grating Spectrometer 
(RGS) exhibits broad [v ~ 1400 km s _1 FWHM) absorption lines from highly ionized elements including 
hydrogen- and helium-like carbon, nitrogen, oxygen, and neon, and several iron L-shell ions (Fexvn - xx). 
Also shown in the spectrum is the first astrophysical detection of a broad absorption feature around A = 16 — 17 
A identified as an unresolved transition array (UTA) of 2p — 3d inner-shell absorption by iron M-shell ions in 
a much cooler medium; a feature that might be misidentified as an Ovu edge when observed with moderate 
resolution spectrometers. No absorption edges are clearly detected in the spectrum. We demonstrate that the 
RGS spectrum of IRAS 13349+2438 exhibits absorption lines from at least two distinct regions, one of which is 
tentatively associated with the medium that produces the optical/UV reddening. 

Key words, atomic processes - line: formation - techniques: spectroscopic - quasars: absorption lines - quasars: 
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1. Introduction 

IRAS 13349+2438 is an archetypal highly-polarized radio- 
quiet quasar at a redshift of z — 0.10764 (Kim et al. 1995| ). 
Since its identification as an infrared-luminous quasar 



(Beichman et al. 1986), this source has been extensively 
studied in the optical, infrared, and X-ray bands. In a 
detailed investigation of the optical and infrared spec- 
tra and polarization measurements, Wills et al. (1992) 
demonstrated that the nuclear spectrum exhibits two dis- 
tinct components; a highly-reddened component and a 
highly-polarized component that suffers much lower ex- 



tinction. Based on these observational facts, Wills (1992) 
constructed a simple and elegant picture of the nuclear 
region of IRAS 13349+2438 in which the direct AGN ra- 
diation is attenutated by a thick dusty torus, while the 



* Based on observations obtained with XMM-Newton, an 
ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and the USA (NASA). 



observed highly-polarized light is produced by scattering 
in an extended bipolar region, either by warm electrons or 
by small dust grains. 

IRAS 13349+2438 was detected in the ROSAT All- 



Sky-Survey (Walter & Fink 1993; Brinkmann & Siebert 



1994J), and has been the target of extensive pointed PSPC 
obser vations with ROSAT (Brandt, Fabia n, fc Pounds 
1996[ ) and with ASCA (Brinkmann et al. |l996[ Brandt 
et al. 1997). The soft X-ray spectrum obtained with the 



PSPC shows a lack of absorption by cold material indi- 
cated by the observed optical reddening, and suggests the 
presence of a warm, dusty medium along the line of sight 
(Brandt, Fabian, & Pounds |1996[ ). 

In a more recent investigation of the complex X-ray 
properties of IRA S 13349+2438, Siebert, Komossa, & 
Brinkmann ( 1999 ) self-consistently accounted for the ef- 
fects of dust embedded in the warm absorbing medium, 
and concluded that single zone models, both with and 
without internal dust, do not provide adequate fits to the 
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combined, ROSAT, ASCA, and optical data sets. In par- 
ticular, they find that a dust-free warm absorber model 
formally gives the best fit to the X-ray data, and con- 
clude that the X-ray absorption and optical reddening 
must arise in spatially distinct regions. However, owing to 
the moderate spectral-resolving-power capabilities of the 
available detectors on ROSAT and ASCA, and the likely 
cross-calibration uncertainties, the precise nature of the 
soft X-ray spectrum has remained controversial. 

In this Letter, we present results from the first high- 
resolution X-ray observation of IRAS 13349+2438 with 
the XMM-Newton Observatory. The spectrum obtained 
with the Reflection Grating Spectrometer (RGS) shows 
a wealth of discrete spectral features, including the first 
astrophysical detection of inner-shell 2p — 3d absorption 
by M-shell iron ions in the form of an unresolved tran- 
sition array (UTA). From a detailed analysis of the rich 
absorption spectrum, we measure the column density and 
velocity field of the line-of-sight material. We show that 
the spectrum contains absorption features from regions 
with two distinct levels of ionization. The column density 
of the lower ionization component is consistent with the 
observed optical reddening, and we tentatively associate 
this component with the dusty torus. 

2. Observation and Data Reduction 

IRAS 13349+2438 was obs erved with the XMM-Newton 
observatory (Jansen et al. 2001 ) on 19 - 20 June, 2000 



during the Performance Verification phase for a total 
exposure time of 42 ks. The data obtained with the 
Refle ction Grating Spectrometer (RGS; den Herder et al. 
2001 ) were filtered through standard event-selection cri- 
teria using the XMM-Newton Science Analysis Software 
(SAS). The source spectrum was extracted using a spa- 
tial filter in dispersion/cross-dispersion coordinates to iso- 
late IRAS 13349+2438 from other possible contaminat- 
ing sources and to reduce contribution from background 
events. Subsequently, the first order events were selected 
by applying a dispersion/pulse-height filter. The back- 
ground spectrum was generated using all of the events 
that lie outside the spatial mask. Wavelengths were then 
assigned to the dispersion coordinates. The current wave- 
length scale is accurate to within ~ 8 mA across the entire 
RGS band of A = 5 - 35 A (E = 0.35 - 2.5 keV). 

The European Photon Imaging Camera (EPIC; Turner 
et al. |200lt Stiider et al. |200ip MOS1 data were also pro- 



cessed with the SAS. The MOS2 detector was operated in 
FAST UNCOMP mode, for which reduction by the SAS is 
not possible as yet. Source events were extracted using a 
circular region of radius 45". A nearby source- free region 
with a radius of 3' was used to assess the background. 

Thre e modes of the Optical Monitor (OM; Mason et 
al. |2000| ) were used during the observation. The V grism 
(or optical grism; exposure time 3000 s), UV grism (expo- 
sure time 1000 s), and the UVW2 filter (~ 1500 - 3000 A; 
effective exposure time 6340 s). There is no significant 
variability in the UVW2 observations. Absolute flux cali- 



brations for the grism data were not available at the time 
of writing. The optical spectrum shows clear evidence of 
H/3 emission but the signal-to-noise is insufficient to de- 
termine the presence of any other emission lines (within 
an observer frame wavelength range of - 3000-6000 A). 
The UV spectrum (coverage ~ 2000 — 3500 A) is too weak 
for spectral extraction. 

3. Results of Spectral Fitting 

3.1. Underlying Continuum Radiation 

To first obtain a rough characterization of the shape of the 
continuum, we use the PN spectrum, which has the high- 
est statistical quality and covers a broad range in energy. 
During the observation, the source exhibited a gradual and 
steady increase in brightness by ~ 30%, but with no noti- 
cable change in the spectral shape. We, therefore, use the 
cumulative spectrum for all of our subsequent analyses. 
The 0.2 — 10 keV spectral region can be well-fitted with a 
phenomcnological model that consists of a powerlaw and 
two black body components absorbed through a Galactic 
colum n density of V| al = 1.1 x 10 20 cm~ 2 (Murphy et al. 



1996). The PN data require two black body components 
with temperatures of kT — 70 eV and — 12 eV. The best- 
fit powerlaw photon index of T — 2.2 is consistent with 
the value implied by the ASCA data (Brinkmann et al. 



1996; Brandt et al. 1997) 



The MOS1 data from 0.3 to 10 keV, excluding the 
0.6—1.2 keV region where the RGS shows complex absorp- 
tion features, were best- fit using a single blackbody plus 
power-law (a % 2 of 1-59 for 199 degrees of freedom), which 
is generally consistent with the PN data. The blackbody 
temperature is kT ~ 110 eV. The best-fit photon index is 
r = 2.1 with a column density of 1.4 x 10 20 cm~ 2 , which 
is only slightly higher than the Galactic value. 

3.2. The Absorption Spectrum of IRAS 13349+2438 

The RGS spectrum shown in Figure [j] exhibits numer- 
ous absorption lines from a wide range in levels of ioniza- 
tion. The most prominent features in the spectrum are K- 
shell absorption lines of H- and He-like carbon, nitrogen, 
oxygen, and neon, and L-shcll lines of Fexvn - xx. The 
spectrum also shows a broad absorption feature between 
A ~ 16 - 17 A (E ~ 730 - 770 eV). The observed location 
of the red "edge" of this feature is at A = 17.10 ± 0.05 A 
in the rest-frame of the quasar, and is close but undoubt- 
edly inconsistent with that of the O vn photoelectric edge 
(A = 16.78 A). The shape of the absorption trough to- 
wards the shorter wavelengths is also incompatible with 
that of a photoelectric edge. We identify this feature as a 
UTA of inner-shell 2p — 3d resonance absorption lines in 
relatively cool, M-shell iron. The shape of this feature is 
strikingly similar to laboratory absorption measurements 
of a heated iron foil (Chenais-Popovics et al. 2000 ), and 
also agrees well with our own calculations as described 
below. 
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Fig. 1. The RGS first order spectrum of IRAS 13349+2438 corrected for cosmological redshift (z = 0.10764). The 
error bars represent la Poisson fluctuations. The wavelength bins are approximately ~ 0.1 A wide. The best-fit model 
spectrum is superimposed in red. Absorption line features labeled in blue are predominantly produced in the low-£ 
component, while those labeled in black originate from the high-£ component (see text for full details). The discrepancy 
at A ~ 29.3 A is probably due to inner-shell absorption by Nv (lithium-like), which is not included in our model. 



We adopt a continuum model similar to that inferred 
from the PN data; i.e., the sum of a powerlaw and a black- 
body component, the latter merely in order to parametrize 
empirically the soft X-ray spectral shape. We also fix the 
powerlaw photon index to T = 2.2 and the column density 
of cold material to the Galactic value. With the continuum 
model defined, we then apply absorption components of El- 
and He-like C, N, O, and Ne, and Fexvn - xxiv. Each ion 
is treated as a separate component in the spectral fit and 
contains all relevant resonance transtions from the ground- 
state and photoelectric edges in the RGS band. The line 
profiles are calculated accounting for thermal and turbu- 
lent velocity broadening. Transition wavelengths and os- 
cillator strengths were calculated with the atomic physics 
package HULLAC (Bar-Shalom et al. 1998), except for 



the wavelengths of the strong Fe L resonance lines where 
we use laboratory measurement values described in Brown 
et al. ( 200C| ). We use photoionization cross sections from 
Verner et al. (1996). For the iron UTA absorption, pho- 



toexcitation cross sections from L-shell to M-shell in Fe v 
- xvi are computed. All of the transitions 2l 8 3l x - 2l 7 



31 



x+l 



1 through 12) are taken into account, of which 



the 2p — 3d excitations are the most important. For the 



atomic structure, the most significant configuration mix- 
ings, which conserve the total angular momentum within 
the n = 3 shell (namely 3p 2 + 3s3d), are included. This ap- 
proximation is expected to be adequate for analysing the 
presently observed unresolved absorption feature. A more 
detailed discus sion of the UTA is presented in Behar, Sako, 
& Kahn (|2000D . 

Each component is then convolved through the in- 
strument line spread function and fit for the ion col- 
umn density simultaneously with the black body con- 
tinuum parameters and the normalization of the pow- 
erlaw component. The best- fit black body temperature 
is kT ~ 100 eV with a flux in the 5 — 35 A range of 
4.0 x 10 -3 photons cm -2 s _1 . The flux in the powerlaw 
component is 2.4 x 10~ 3 photons cm~ 2 s _1 in the same 
wavelength range. 

The observed widths of the absorption lines, as shown 
below, are much larger than both the instrument line 
spread function and the thermal broadening of a gas with 
kT ~ 10 eV, which is the expected temperature for a 
photoionized plasma at this level of ionization. With the 
current statistical quality of the RGS data, however, we 
are not able to constrain the turbulent velocities v tur h of 
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the individual ion components. We, therefore, assume a 
uniform mean turbulent velocity field, keeping in mind 
that each ion can, in principle, exist in regions of differ- 
ent turbulent velocities. The derived ion column densities, 
therefore, may be uncertain to some degree, as quantified 
in the following section. 

We obtain a statistically acceptable fit to the RGS data 
with Xr = 1-07 for 532 degrees of freedom. The continuum 
parameters inferred from the RGS data are also consis- 
tent with those derived from PN and MOS. The intrinsic 
isotropic luminosities in the 0.3-2 keV and 2-10 keV 
regions are Lx ~ 2 x 10 44 erg s _1 and ~ 5 x 10 43 erg s _1 
(assuming H = 65 km s _1 Mpc -1 and qo = 0.5), respec- 
tively, which are lower than both the ROSAT PSPC and 
A SCA values by a factor of ~ 5. The measured ion column 
densities are listed in Table |[ To illustrate the statistical 
significance of the various components, we also list the 
changes in \ 2 when each of the components is removed 
from the model. The best-fit FWHM turbulent velocity is 
^turb — 1430^280 s_1 > which is much larger than the 
thermal velocity of a photoionized medium. We also find 
weak evidence of an average bulk outflow velocity shift 
with Ughift = +200^180 km s _1 , where positive velocity 
denotes a blueshift. 



Table 1. Measured Ion Column Densities 



Ion 


N t ( 


cm 


2 ) 


a 


A X 2b 


Component 


C V 


(6.3 


+ 5.2n 
-4.8/ 


X 


10 ib 


4.8 


low 


C VI 


(6.0 


+6.4\ 
-3.1/ 


X 


10 16 


23.7 


low 


N VI 


(2.4 


+1.4\ 
-1.0/ 


X 


10 16 


29.1 


low 


N VII 


(1.3 


+0.8 \ 
-0.4) 


X 


10 17 


109.5 


low/high 


O VII 


(3.7 


+2.7 \ 
-1.8/ 


X 


10 16 


24.2 


low 


O VIII 


(9.5 


+8.7\ 
-4.7/ 


X 


10 16 


24.7 


high 


Ne IX 


(1.2 


+0.6 \ 
-0.5/ 


X 


10 18 


9.9 


high 


Ne X 


(4.9 


+ 10 \ 
— 3.2/ 


X 


10 17 


25.9 


high 


Fe XVII 


(1.7 


+ 1.6 \ 
-1.2/ 


X 


10 17 


23.0 


high 


Fe XVIII 


(6.3 


+2.4\ 
-1.9/ 


X 


10 17 


115.0 


high 


Fe XIX 


(9.7 


+4.8% 
-3.9) 


X 


10 17 


166.7 


high 


Fe XX 


(3.0 


+2.8 \ 
-2.0) 


X 


10 17 


39.5 


high 


Fe VII 


(1.5 


+1.5\ 
-1.3/ 


X 


10 16 


9.6 


low 


Fe VIII 


(4.6 


+1.5\ 
-1.3^ 


X 


10 16 


91.8 


low 


Fe IX 


(8.8 


+ 12 \ 
— 8.7/ 


X 


10 15 


29.3 


low 


Fe X 


(2.4 


+ 1-3-1 
-1.1/ 


X 


10 16 


72.4 


low 


Fe XI 


(1.9 


-0.9) 


X 


10 16 


65.3 


low 


Fe XII 


(6.4 


+ 10 \ 
-4.4 ! 


X 


10 16 


25.6 


low 



a The turbulent velocities for all of the ions are fixed to the 
same value. We find a best-fit with Uturb = 1430^280 km s _1 
FWHM. 

b The increase in x 2 when the ion component is excluded from 
the best-fit model. 

c The dominant ^-component responsible for producing the 
absorption feature (see Figure |l| and text for details). 



4. Implications of the Results of Spectral Fitting 

For the measured ion column densities listed in Table [l], 
many of the strong absorption lines of neon and iron L ions 
are in the logarithmic region of the curves of growth, which 
indicates that the derived column densities are highly cou- 
pled with the assumed turbulent velocity. Those of K-shell 
carbon, nitrogen, and oxygen, and iron M-shell ions, how- 
ever, are in the linear regime, and the derived column 
densities are not very sensitive to the turbulent velocity. 
An important point to note is that the observed veloc- 
ity widths may as well be due to a superposition of mul- 
tiple, discrete absorption components each of which are 
optically thin, and are unresolved with the RGS. Such 
features have been observed in UV absorption line s pec- 
tra of many Seyfert 1 galaxies (Crenshaw et al. 1999 and 
references therein), some of which show as many as 7 dis- 
tinct, kinematic components (e.g., Mrk 509; Kriss et al. 
[2000D . If this is the case, the measured column densities 
from the neon and iron L lines, which are produced in the 
logarithmic region of the curve of growth, may be under- 
estimated relative to the true value. The column densities 
of the He-like species, as well as those of the UTA, are not 
affected by this, since the total observed values are still in 
the optically thin regime. 

As shown in Table the detections of absorption lines 
from Fevn - XII and Fexvn - xx are highly significant. 
The column densities of the intermediate charge states 
of Fexiii - XVI, however, are consistent with zero, with 
an upper limit of Ni ~ few x 10 15 cm -2 for each of 
these charge states. This indicates that the line-of-sight 
material consists of either a multi-phase gas in a sin- 
gle medium, or two (or more) spatially distinct regions. 
Motivated by this observational fact, we refit the spec- 
trum using the same model, except, we assume that the 
line-of-sight material consists of two discrete velocity com- 
ponents; (1) the low-ionization-parameter component in- 
cluding C vi, C v, N vi, O vii, and the M-shell Fe, and (2) 
the high-ionization-parameter component including Nvn, 
O VIII, Neix, Nex, and L-shell Fe. Contrary to our previ- 
ous fit where the bulk velocities of all the ions were fixed 
relative to one another, we find that the low-ionization- 
parameter component is significantly blue-shifted with 
Vshift.iow = +420j^gQ km s" 1 , while the bulk velocity shift 
in the high-ionization component is consistent with zero 
(w s hift,high = —201.330 km s_1 )- This implies that the low- 
ionization gas is being accelerated substantially compared 
to the high-ionization component, as would be expected 
in a radiatively-driven outflow (Arav & Li 1994). The de- 



rived turbulent velocities of both components remain con- 
sistent with that of our previous fit (i>t ur b ~ 1500 km s _1 
FWHM). 

From the observed distribution of charge states of M- 
shell iron, the average ionization parameter, £ = L/nr 2 
(erg cm s _1 ), of the absorbing gas is estimated to be 
log£ ~ based on a calculation with the photoionization 
code XSTAR (Kallman & Krolik |l995[ ) using the inferred 
continuum shape for the ionizing spectrum. The width of 
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the distribution in ionization parameter is no larger than 
A log £ ~ 1 . The measured iron ion column densities sug- 
gest that the corresponding equivalent hydrogen column 
density is JVh ~ (1 — 3) x 10 21 cm~ 2 assuming a solar iron 
abundance. This low-£ gas accounts for most of the car- 
bon and He-like nitrogen and oxygen absorption lines as 
well, however, with some indication of lower carbon and 
nitrogen abundances by a factor of ~ 2 and oxygen by a 
factor of ~ 3 relative to the solar iron abundance. The lack 
of absorption from Fexni - xvi, however, indicates that 
substantial amounts of material with ionization parame- 
ters in the range 1 55 log£ Sa 2 are not present along the 
line-of-sight. Whether this is related to the global struc- 
ture of the circumnuclear medium or a mere coincidence 
from a superposition of physically distinct regions is not 
known. 

The absorption lines from H-like nitrogen and oxygen, 
H- and He-like neon, and L-shell iron are produced in a 
medium with 2.0 ^ log£ ^2.5 and an equivalent hydro- 
gen column density of ATjj ~ (1 — 4) x 10 22 cm~ 2 , again, 
assuming a solar iron abundance. The inferred abundances 
of nitrogen, oxygen, and neon relative to that of iron are 
consistent with solar values, but are not well-constrained. 

For a normal dust-to-gas ratio, the observed redden- 
ing of E (B - V) = 0.3 in IRAS 13349+2438 (Wills et 
al. 1992 ) corresponds to a hydrogen c olumn density of 
N H = 1.7x 10 21 cur 2 (Burstein & Heiles [L978|). This value 



is significantly lower than the total amount of X-ray ab- 
sorbing material (low-£ + high-£) observed in the present 
X-ray spectrum. Coincidentally, however, the derived col- 
umn density of the low-£ component is very close to that 
of the optical reddening, although, we cannot conclusively 
associate the low-£ X-ray absorber with the dusty torus. 
The column density of the high-£ component, on the other 
hand, is a factor of ~ 10 higher. 

An interesting point to note is that the 1 55 log£ ^ 2 
region is not thermally unstable, based on XSTAR cal- 
culations described above. On the other hand, the high-£ 
region (2.0 ^ log£ Si 2.5) that we observe in the spec- 
trum is thermally unstable. However, complications such 
as non-solar metal abundances and/or inaccurate ioniza- 
tion and recombination rates may alter the shape of the 
thermal stability curve significantly, and, hence, the tem- 
perature ranges of the unstable regions (Hess, Kahn, & 
Paerels |l997| ; Savin et al. |1999p . 

For the observed X-ray luminosity of Lx ~ 2.5 x 
10 44 erg s _1 , the high-£ component with log£ ^ 2.0 pro- 
vides the constraint, n e r 2 < 2.5 x 10 42 cm -1 . The mea- 
sured column density through this medium is Nh ~ 
n e Ar il 10 22 cm -2 , where Ar is the radial thickness. 
Assuming that Ar/r < 1, these constraints combined pro- 
vide an upper limit in the location of the high-£ gas of 
^high-C < 2.5 x 10 20 cm ~ 80 pc, which is representa- 
tive of a typical narrow-line region. The correponding gas 
density in this region is n e <; 40 cm -3 with an estimated 
thickness of Ar ~ 80 pc ~ rhi g h-£- A similar calculation 
for the low-£ component does not provide a useful con- 
straint with < 10 23 cm ~ 3 kpc, and, therefore, the 



location of the low-£ component is not well-determined 
compared to that of the high-i; component. If the low-£ 
component lies beyond the high-£ component relative to 
the central continuum source, the difference in the mea- 
sured bulk velocity shifts of the two components might 
indicate that the high-£ gas is the base of an accelerat- 
ing outflow. If, on the other hand, the low-£ absorber is 
indeed spatially coincident with the torus, in which case 
1ow-£ is likely to be approximately the location of the 
broad- line region {r\ ow -£ ~ ?*blr ~ 1 P c ) ; the medium is 
decelerating as a function of radius. Such a behavior has 
been observed in the UV spectrum of the Seyfert 1 galaxy 
NGC 4151 (Crenshaw et al. |200C| ). 

As briefly mentioned earlier, the source during the 
present XMM-Newton observation was in an unusually low 
state. However, since the estimated location and density 
of the absorbing medium are such that the gas does not 
respond immediately to the observed continuum radiation 
(i.e., low density gas at large distances), the effect on the 
physical state of the absorbing medium by a variable illu- 
minating source is not clear. It will be useful to re-observe 
IRAS 13349+2438 during a substantially brighter state to 
see whether the spectrum exhibits any dramatic changes 
in the absorption structures, and specifically if the oxygen 
absorption edges detected i n the ROSAT and A SCA data 
(Brandt, Fabian, & Pounds |l996| ; Brandt et al. |l997|) are 
really present during a different state. A detailed compar- 
ison of the UV and X-ray absorption spectra will also be 
interesting, particularly for identifying discrete kinematic 
components, as well as for constraining the global charac- 
teristics and dynamics of the absorbing medium. 

5. Comparisons with Other AGNs 

The absorption spectrum of IRAS 13349+2438 is qualita- 
tively similar to those obtained with the Chandra trans- 
mission grating observations of the Seyfert 1 galaxies NGC 
5548 and NGC 3783, which show narrow absorption lines 
blue-shifted by several hundred km s^ 1 . (Kaastra et al. 
|20U0| ; Kaspi et al. |2000| ). The derived ion column densities 
in these sources, as well as in IRAS 13349+2438, are in 
the range N% ~ 10 16 — 10 17 cm -2 , and are not high enough 
to produce observable absorption edges. 

Conceptually, the low-ionization component observed 
in IRAS 13349+2438 is similar to the "lukewarm ab- 
sorber" that explains both the observed optical and X- 
ray attenuation in NGC 3227 (Kraemer et al. 2000). 



The spectroscopic signatures, however, are very different. 
In particular, for the column densities derived from the 
IRAS 13349+2438 data, the absorption features are dom- 
inated by discrete resonance line transitions, mainly in 
He-like ions and M-shell iron, and not by ph otoelectric 
edges as in the model of Kraemer et al. ( 2000 ). 

As demonstrated in our detailed spectral analysis of 
IRAS 13349+2438, the UTA feature is potentially a pow- 
erful diagnostic tool for probing cool absorbing mate- 
rial using high-resolution X-ray observations. If the low- 
ionization component is indeed associated with the dusty 
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torus as the derived column density suggests, this feature 
should be detectable in other AGNs where the line-of-sight 
is partially obscured by the torus. 
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